Size-selective nanoparticle growth on few-layer graphene films 
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ABSTRACT 

We observe that gold atoms deposited by physical vapor deposition onto few layer 
graphenes condense upon annealing to form nanoparticles with an average diameter that is 
determined by the graphene film thickness. The data are well described by a theoretical model in 
which the electrostatic interactions arising from charge transfer between the graphene and the 
gold particle limit the size of the growing nanoparticles. The model predicts a nanoparticle size 

distribution characterized by a mean diameter D that follows a D cc m 1/3 scaling law where m is 
the number of carbon layers in the few layer graphene film. 
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Graphene, a honeycomb sheet of sp -bonded carbon atoms, has attracted immense 
attention since the first demonstration of reliable methods to rapidly isolate and identify it. 1 More 
recently it has become clear that graphene multilayers or "few layer graphene" (FLG) films 
actually represent a family of materials with properties that vary systematically with the number 
of carbon layers. As an example, prior theoretical work considered the electronic structure of 
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FLG and the effect of a substrate on the spectrum. " More recently, electrostatic force 
microscopy was used to measure very significant variations of the surface potential with 
thickness for FLG films on an oxidized silicon substrate, which were found to be in excellent 
agreement with a nonlinear Thomas-Fermi theory for the FLG charge carriers. Here we report on 
another important effect of FLG film thickness, namely its impact on the diameter of nanoscale 
metal particles formed on the FLG surface. 

Thick graphite samples are widely used as substrates for the preparation of free-standing 
metal nanoparticles due to their chemical inertness, low binding energies for adsorbed atoms, 
low energy barrier for lateral diffusion, and suitability for various electron microscopies, in 
particular for scanning tunneling microscopy. ' Here we study the formation of nanoparticles 
on FLG films, and observe that the size of the growing nanoparticles is controlled by film 
thickness. We find that mean particle diameter D increases with the layer count of the graphene 
m, as does the variance of the distribution of particle radii. The observations are compared with 
the predictions of a theoretical model where long-range electrostatic interactions resulting from 
charge exchange between graphene and the nanoparticle limit the growth of large nanoparticles. 
The nanoparticle diameter is determined by a competition between its electrostatic energy and 
the surface tension of the particle. The model predicts that D increases as floe m 1/3 ; while the 
variance of the radius distribution grows as^~«. The power law dependence of D occurs as 
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a consequence of the unconventional charge screening in the graphene film. 

We find very good agreement between the data and the predictions of the model. 
Specifically, we fit the average diameter data to a power law D= Am v ; as suggested by the 
model- The best- fit values for the exponent and prefactor are v= 0.33 ±0.06 anc j 

A — 6.46 ± 0.68 nm ? respectively, which compare well with the theoretical predictions ^=1/3 
and ^ = 5-9 nm. Good agreement is also found for the variance data (see below). 

Graphene samples were obtained by micromechanical cleavage of Kish graphite and 
placed on 300nm oxidized substrates. 1 The number of graphene layers was determined by cross- 
corroboration between interference induced color shift in optical microscope and atomic force 
microscope (AFM), and further confirmed by their Raman spectra 11 ' 12 (see Fig. SI of the 
Supplemental Information). Sample cleanliness was found to be an important factor in obtaining 
reproducible results, so the as-prepared graphene samples were annealed at 400°C in forming gas 
(10% hydrogen, 90% ultrahigh purity argon gas) for 5 hours to remove adsorbates and 
contamination. Thermal evaporation of gold was performed using an alumina-coated tungsten 
boat at vacuum pressure of 5x10" Torr, while the deposition flux (0.05 nm/s) and final thickness 
(0.3 nm) were monitored with a quartz crystal microbalance. The sample was then loaded in a 
quartz tube chamber and annealed at 400°C for 3 hours in a flow of forming gas. This allowed 
nanoparticles to form and reach an equilibrium configuration, as evidenced by the fact that no 
further particle morphology change was observed when such a sample was annealed for 4 
additional hours. 

Figures 1(a) and 1(b) are high-resolution scanning electron micrographs of identically 
prepared samples of Au nanoparticles that differ only in the thickness of the FLG film (~ 700 
layers for Fig la and a single layer for Fig lb; see Supplemental Information). Striking 

3 



differences are observed in the size and shape of the gold nanoparticles between thick and thin 
films. 

For thick FLG films on oxidized silicon substrates, we find a nanoparticle morphology 
(Fig la) identical to that reported for gold nanoparticles formed on bulk HOPG after vacuum 
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annealing, indicating that the presence of the silicon substrate is negligible in this configuration. 
There are a comparatively small number of large nanoparticles (diameter ~ 100-300 nm) whose 
facets are oriented along directions rotated by multiples of 60 degrees, as indicated in Fig la and 
analyzed in the inset. This is similar to what is observed in the epitaxial growth of micron-sized 
Au particles on graphite 14 where the {111} face of the metal grows on the {0001} face of 
graphite. In this case, particle growth is governed by the well-studied Ostwald ripening process, 
where large particles grow at the expense of the smaller particles, because of the latter 's 
proportionately higher surface energy. 15 

For Au deposited onto single- and few-layer graphene films (Fig. lb), however, the 
nanoparticle morphology after annealing is radically different. The FLG film surface is now 
decorated by a large number of comparatively small nanoparticles (diameter 5-20 nm) with no 
evidence of faceting. This configuration represents the equilibrium state of the system since no 
significant further size changes are observed after several additional hours of annealing. 

We find that the average nanoparticle diameter, measured by HRSEM, and the statistical 
variation in this quantity are strongly dependent on the thickness of the FLG film, with the 
smallest particles and narrowest distribution observed on single layer graphene (see, for example, 
Fig. 2a which shows particles grown on FLG films of 1, 2, and 3 layers thickness). These 
quantities are plotted in Fig. 2b and its inset, along with power law fits that are motivated by an 
analysis of the particle free energy as discussed below. Additional detail and AFM and SEM 
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images are provided in Figs. S5-12 in the Supplemental Information. We also find that the 
average diameter-to-height ratio for nanoparticles grown on FLG films is approximately 5, so 
that the particles grow as relatively "flat" islands. This is inferred from Fig. 3, which shows the 
height (measured by AFM) versus diameter for nanoparticles grown on a FLG film. For this plot, 
the particle height is measured by AFM (accuracy of 0.1 nm), while the lateral dimension of the 
same particle is measured by HRSEM in order to avoid tip convolution broadening of 5-10 nm 
associated with the AFM tips used (see Fig. S4 in Supplemental Information for details). Thus, 
nanoparticles grown on FLG films of fewer than 20 layers grow as relatively flat cylinders, with 
an average diameter controlled by the FLG thickness. 

Our observation that the equilibrium configuration of metal nanoparticles formed on thin 
FLG films differs strongly from that of particles grown on thick graphene films suggests that the 
interactions that control particle growth differ sharply as well. The large particles on thick 
graphene apparently grow by the Ostwald ripening process that is determined by the interplay of 
bulk cohesive energy and surface interactions. The different morphology of nanoparticles on 
FLG indicates the impact of a new long-range interaction whose strength grows with decreasing 
FLG thickness. Prior work ' suggests the likelihood of an electrostatic interaction whose 
strength is thickness dependent. Indeed, we find that the data are well described by a model 
where electrostatics play a critical role in the size selection, and so we examine this interaction in 
detail before constructing a model for the size-dependent total energy of the particle. 
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Consider the effect of the work function mismatch A<P between graphene (-4.66 eV ) 
and gold (5.1-5.47 eV). When the gold nanoparticle is adsorbed on graphene there is a transfer of 
electrons from the graphene to the nanoparticles, leaving the particle negatively charged with an 
interfacial charge density <r . Each donated charge produces a dipole of strength p = es where s 
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parameterizes the effective width of the dipole layer. This width is determined by the physical 
penetration depth of the interfacial charge density on the gold particle and the depth of the 
screening charge in the graphene film. Since few layer graphene films screen poorly, s is 
effectively limited by the penetration depth in the graphene. An areal dipole density t = A 0/4 

n e s is then required to equilibrate the Fermi energies of the metal nanoparticle and the 
graphene. 

The interfacial dipoles are all equal and oriented normal to the surface with the positive 
end outward, leading to a (positive) electrostatic energy for the droplet 

AO 2 rd 2 rd 2 r' 
327T e J \r — r \ 

where the integral is taken over the cross sectional area of the droplet. This integral is formally 
divergent because of the small distance behavior of the integrand. The physical energy is finite 
since the integral is cutoff at small length scales by the spatial extent of the dipole. Below this 
scale (1) reverts to a nonsingular near field form that makes a contribution to the bulk energy, 
proportional to the droplet volume. Thus, by integrating (1) for r >d, where d is the thickness 
of the graphene film, we obtain the long-range contribution to the energy that shows a new 
scaling with droplet size. 

The interaction energy can be equivalently expressed, using a Fourier transform, as 
U d = T^\d 2 qV(q)\S(qf (2) 

f R 

where S(q) = 2n) Q rJ {qr) dr is the form factor for a circular island with radius R. By cutting off 

the interaction range at the graphene film thickness we find that V(q — » 0) = 1 / d and thus the 
long range electrostatic interaction energy becomes 
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The areal density r appearing in this expression depends on the Fermi energy mismatch A <X>, and 
it can be deduced from a Thomas Fermi model for interlayer screening developed by us 

previously in Ref. 8 which gives r = CA<1> 3 2 with C = V3/8 {l/ ftv F 4e^cy where c = 0.34 nm 

is the interlayer spacing in graphite. Thus, eliminating the areal dipole density, the electrostatic 
energy reads 

U d = 4-z^Tr^ 1!2 = L R A (4) 
" 8 3/2 hv F de 3 ^ m 

where m = dl c is the layer count of the graphene. This is a useful result, since it demonstrates 

that the long range electrostatic energy scales as the square of the droplet area, inversely with the 

film thickness and vanishes if there is no charge transfer because of matched work functions. 

Assuming a work function difference AcD = 0.5 eV, we find a value for Y = 0.61 eV/nm 4 . 

This electrostatic interaction can now be included to complete a model for the total 

energy of a condensed droplet. This energy has contributions from the particle surface area A 

and its volume V as well as the electrostatic contribution 

U=yA + uV+-R 4 (5) 
m 

For a cylindrical particle with radius R and height h we have A = 2nR 2 + 2nRh , V = nR 2 h, y 

represents the surface tension (energy/area; from the literature 18 , we find y = 9.6 eV/nm 2 ) and u 
(<0) is the bulk cohesive energy per unit volume. In principle, one needs to treat separately the 
energies of the free and contacted (bottom) surfaces of the droplet, but this detail is unnecessary 
in the following scaling analysis. We note that in the absence of the electrostatic term, the model 
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exhibits the familiar critical droplet phenomenon at a critical radius R c = yhl (\u\h- 2y) below 

which the particles collapse and above which the particles grow without bound. The long-range 
electrostatic term then has the effect of limiting the growth of large droplets. It is useful to 
express the result by computing the chemical potential 



where n is the atomic number density (for gold, n = 59 / nm ), and u 1 = u + T I h gives the 
(constant) volume dependent term in the droplet energy. 

Figure 4(a) illustrates a plot of the phenomenological free energy as a function of disk 
diameter for various thickness of graphene, as obtained from Eq. 6. The total droplet energy 
increases linearly for small-diameter particles because its free energy is dominated by the line 
tension on the boundary. For intermediate particle diameters, the negative volume energy makes 
the overall energy negative. Electrostatic energy dominates the large R behavior of the energy 
and prevents a runaway growth of the droplet. Figure 4(b) illustrates chemical potential of 
nanoparticles grown on different number of graphene layers, calculated using Eq 6. The 
chemical potential diverges <x 1/R for small R from the surface tension term and it increases <x 
R at large R because of the dipolar interactions. The chemical potential curve sharply rises for 
small diameters, with a gradual increase for large diameters. The curve has a broader minimum 
for thicker graphene, which predicts a broader diameter distribution for such flakes. 

In equilibrium the system forms droplets at the average radius that minimizes the 
chemical potential /j(R) . Since the potential of Fig. 4(a) is concave downward, in practice this 
implies a microphase separation of the adsorbed gold into a low density (uncondensed) phase 
and droplet (condensed) phase at an average radius corresponding to this minimum. From Eq. (6) 




(6) 
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we can predict that for an ensemble of droplets in thermal equilibrium, the average diameter 
follows: 



D = 2R = 



= 5.9m m nm (7) 



V r ) 

The curvature of ju(R) around its minimum value is /u" = 12r / nhmn . Thus for an equilibrium 
ensemble of droplets that are formed at the annealing temperature T , one predicts a distribution 
of droplet radii with a variance that satisfies: 
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(8) 



ol = R - R = = 1.6m nm 2 

R 12T 

where { ) denotes a thermal average. To determine the numerical prefactors in Eq. 7-8, we use 
the values of y , T, and n given previously, along with a value of h = 2 nm taken from the data 
in Fig. 3. 

These theoretical predictions motivate the curve fits to the data presented in Fig. 2(b) and 
its inset. The particle sizes and standard deviations increase with graphene thickness. Fitting the 
average diameter data to a power-law yields an exponent of 0.331 ± 0.061 and a prefactor of 6.46 
± 0.68 nm, in good agreement with the theoretical predictions of 1/3 and 5.9 nm, respectively. 
The measured variance of the diameter distribution is presented in Fig. 2(b). The variance is 
found to increase linearly with the number of graphene layers, in agreement with the theoretical 
prediction. The best fit slope is 0.71 nm 2 , about a factor of 2 below the prediction of the theory, 
and there is a constant offset of 3.67 nm that we ascribe to uncertainty associated with the 
HRSEM measurement of the particle diameter. 

In conclusion, we report a new method for synthesizing metal nanoparticles of selected 
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size. This involves physical vapor deposition of metal onto graphene on SiCVSi wafers, followed 
by aging toward equilibrium in forming gas. Our modeling suggests that this method utilizes the 
differential electrostatic energy of graphene due to its strong thickness-dependent interlayer 
charge screening arising from the relativistic low energy carrier characteristic. The repulsion of 
dipoles induced in the nanoparticles limits the size of the formed nanoparticles. Measurements of 
the average nanoparticle diameter and the variation of this quantity are in very good agreement 
with theoretical predictions that the optimum nanoparticle size follows a simple 1/3 power law 
scaling rule with graphene thickness, and a linear increase in the variance of the size distribution. 
This method could be adapted for the growth of nanoparticles with controlled dimensions 
suitable for application in electronics, optical and magnetic devices, catalysts, and 
nanometrology. 
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FIG. 1. HRSEM images of Au nanoparticles grown on graphite with thickness of 1.0 urn and a 
single layer graphene flake. Inset in Fig. 1(a) indicates the histogram of nanoparticle edge 
orientations. 
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FIG. 2. (a) HRSEM images of Au nanoparticles prepared on few layer graphene on SiC^/Si 
substrate. The numbers and arrow indicate the number of graphene layers. Inset in (a) is the 
AFM image of the same region before nanoparticle growth, (b) Average diameter of Au 
nanoparticles as a function of the number of graphene layers. The red line is a power law fit 
suggested by theory. Inset in (b) is the corresponding dependence of the variance of the diameter 
distribution; the curve shows the theoretically predicated linear dependence. The data point for 9 
layers is shown in the graph but omitted from the fit. 
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FIG. 3. Plot of the height of individual Au nanoparticles as a function of their lateral diameters. 
The inset shows a schematic of the sample. 
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FIG. 4. (a) Plot of the phenomenological free energy as a function of disk diameter for various 
graphene thicknesses. The inset is an expanded view of the region in the dashed rectangle, (b) A 
plot of the chemical potential as a function of the disk diameter for graphene of various layer 
numbers. The inset is a schematic view of the system studied. 
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